Statistical characterization of the forces on spheres in an upflow of air 
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The dynamics of a sphere fluidized in a nearly-levitating upflow of air were previously found to 
be identical to those of a Brownian particle in a two-dimensional harmonic trap, consistent with a 
Langevin equation [Ojha et al, Nature 427, 521 (2004)]. The random forcing, the drag, and the 
trapping potential represent different aspects of the interaction of the sphere with the air flow. In 
this paper we vary the experimental conditions for a single sphere, and report on how the force 
terms in the Langevin equation scale with air flow speed, sphere radius, sphere density, and system 
size. We also report on the effective interaction potential between two spheres in an upflow of air. 

PACS numbers: 05.10.Gg, 47.27.Sd, 47.55.Kf 
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I. INTRODUCTION 

One of the great challenges in physics today is to 
understand the dynamics of driven nonequilibrium sys- 
tems jlj. This is particularly important in soft-matter 
physics, because the materials often have a delicate meso- 
scopic structure that is easily perturbed far from equi- 
librium. There, an understanding of the microscopic 
dynamics is crucial for a fundamental understanding of 
macroscopic behavior. Granular materials are an excel- 
lent example of this point 0, 01 ■ When subjected to 
strong driving forces, granular systems exhibit gas- or 
liquid-like behavior at the macroscopic scale and strong 
velocity fluctuations and collisions at the grain-scale. 
The microscopic fluctuations are created by the act of 
flowing, and, at the same time, are responsible for the 
dissipation that limits the rate of flow. The difficulty 
of treating the fluctuations is one reason why granular 
mechanics remains a forefront research topic, and why 
engineering systems are alarmingly prone to failure. 

One way to characterize the microscopic dynamics in a 
granular gas or liquid is by the distribution of speed fluc- 
tuations. This has a long history, and is associated with 
attempts to develop a system of partial differential equa- 
tions describing granular hydrodynamics [1, IM Ifl 0> la 13 • 
The average kinetic energy associated with speed fluctu- 
ations has come to be known as the "granular tempera- 
ture" , in loose analogy with kinetic theory of gases. An 
interesting line of research has been to explore the ex- 
tent to which this analogy holds, i.e. the extent to which 
statistical mechanical concepts for true thermal systems 
can be used to describe granular fluctuations. For di- 
lute or two-dimensional systems it is relatively straight- 
forward to track grain motion by video techniques. Ex- 
perimentalists have thus studied whether or not speed 
distributions are Gaussian , an d whether or not equipar- 
tion is obeyed HHIl [H llllll . Recently we did the 
same for a very dilute system, consisting of only a sin- 
gle grain, driven by a steady upflow of gas 15|. Part 
of our motivation was to isolate the role of gas-mediated 
interactions from coUisional and cohesive interactions in 
bulk gas-fluidized beds, which is a topic of long-standing 
importance jla| . By measuring the time-dependent dy- 



namics, as well as the usual speed distribution, and by 
performing auxiliary mechanical measurements, we were 
able to demonstrate that the motion of the sphere is iden- 
tical to that of a Brownian particle in a two-dimensional 
harmonic trap. For such a system the thermal analogy is 
perfect. 

In this paper we exploit the thermal analogy to de- 
duce quantitative information about interactions in gas- 
fluidized systems. Now that the tools of statistical me- 
chanics are at our disposal, we may deduce the salient 
features of the forces acting on a sphere from measure- 
ments of position and speed statistics. Besides provid- 
ing additional data and a more detailed description than 
in Ref. |l5|, this follows through on our original moti- 
vation to study the fundamental forces at play in gas- 
fluidized beds. Our statistical mechanical approach is 
completely orthogonal to traditional wind-tunnel mea- 
surements )l7l |. and provides a clean decomposition of 
gas-mediated interactions into three distinct contribu- 
tions. We begin with a discussion of statistical mechan- 
ics and the Langevin equation of motion, both to review 
prior findings and to define notation for use here. After 
describing our experimental apparatus, we then present 
data pertaining to the effective temperature and its scal- 
ing with system parameters, all for a one-sphere system. 
Lastly, we turn to interactions of a sphere with both the 
container boundary as well as with a second sphere. 



II. LANGEVIN EQUATION 

The particles of interest are spheres of mass m, diame- 
ter D = 2Rd, and moment of inertia /. They roll without 
sliding, so their kinetic energy \s K = ^{m+I/Rj^)v'^. In 
order to characterize the motion entirely in terms of posi- 
tion, velocity, and acceleration vectors, {r(t), v(t), a(i)}, 
we define an effective inertial mass and density as m^ — 
m + I/Rd' andpe = me/[(47r/3)i?d'^], respectively. As 
shown in Ref. |l5J | , the equation of translational motion 
of the rolling gas-fiuidized sphere is 



mea(i) = — Vy(r)— m^ 



T{t-t')^r{t')dt'+^Fr{t). (1) 



This is recognized as Newton's Second Law, where the 
right hand side is the sum of forces acting on the sphere. 
The first term is the gradient of an effective potential; 
for a harmonic spring this force is —Kr{t). The second 
term represents the drag force, where r(i) is the memory 
kernal. In Ref. jlfi| it was shown to be exponential. 



r(i) = ro7oCxp(-7ot). 



(2) 



Thus 1/70 is a time scale representing the duration of the 
memory; l/Fo is a time scale such that the drag force has 
a typical value of — TOeFoV. The final term in Eq. Ip^l is a 
time- varying random force Fr{t). As shown in Ref. |15J . 
the components of Fj.(i) have Gaussian distributions and 
exponential temporal autocorrelations. In particular, it 
was demonstrated that the random and drag forces are 
related according to the Fluctuation-Dissipation Relation 
(FDR) [13: 



{Frit') ■ F,.(i)) = 2mekTr{t ~ t'), 



(3) 



where kT = me{v'^)/2 is the effective (granular) temper- 
ature. Note that the two momentum degrees of freedom 
each have kT/2 of energy, consistent with the equipar- 
tition theorem. Satisfaction of the FDR means that 
the particle dynamics are identical to those of a ther- 
mal Brownian particle; therefore, Eq. (Q is truly the 
Langevin Equation. Even though the rolling sphere is 
a driven far-from-equilibrium system, statistical mechan- 
ics holds unchanged except that the value of the effective 
temperature is not the thermodynamic temperature of 
the apparatus. 



III. EXPERIMENTAL DETAILS 

Our methods are identical to those first reported in 
Ref. 15j. The fluidization apparatus is built around a 
12-inch diameter brass sieve, with 300 /im wire mesh 
spacing and with 4 inch high side wall. The full sieve 
is usually used, but occasionally a cylindrical insert is 
placed concentrically in order to vary the radius Rceii 
and/or the wall height. The wire mesh is fiat and level, 
and is very fine compared to the sphere size. The sieve is 
mounted atop a 20 inch x 20 inch x 4 foot tall windbox 
consisting of two nearly cubical chambers separated by a 
perforated metal sheet. In some of the runs, a 1/2-inch 
thick foam air filter is sandwiched between a second per- 
forated metal sheet. Air from a blower is introduced to 
the lower chamber through a fiexible cloth sleeve. The 
fiow rate is controlled by a variac. The geometry of the 
windbox is designed to achieve a uniform upwards air 
flow across the whole area of the sieve. This is verified 
and monitored with a hotwire anemometer. 

The sphere position is measured from digital images 
acquired at a rate of 120 frames per second. The cam- 
era has a resolution of 640 x 480 pixels, and is mounted 
about 1 meter directly above the sieve via a scaffolding 
attached to the windbox. Two 18-inch fluorescent lights 



are mounted just below the camera, such that the illu- 
mination is uniform and the thresholded image of the 
sphere appears as a white disk in a black background. In 
order to achieve very long run times, using an ordinary 
personal computer, we developed custom video compres- 
sion and particle tracking algorithms that permit real- 
time analysis without the need for writing prohibitively 
large data-sets to hard-drive. At heart is a run-length 
encoding scheme: for each row, it's enough to note the 
starting pixel and the segment length. Since black pixels 
have zero intensity, the sphere location is then computed 
as the center of brightness of the entire thresholded im- 
age. 

The sphere velocity and acceleration are found by post- 
processing position vs time data. Specifically, we fit a 
third-order polynomial to data within a window of ±4 
points. Gaussian weighting that is nearly zero at the 
edges is used to ensure continuity of the derivatives. This 
process also refines the position measurement. In the end 
we achieve a resolution of ±0.05 mm, which corresponds 
to about 0.1% of the sphere diameter and about 0.08 
pixels. 

The specific spheres studies are listed in Table 13 For 
each, the allowed air speeds u are bounded by 200 and 
500 cm/s depending on the sphere. The range is lim- 
ited because at lower air speeds, the sphere occasionally 
rolls along its seem or along the weave of the wire mesh. 
At higher air speeds, the sphere occasionally scoots or 
loses contact with the sieve. In all cases, the air speed 
is less than the terminal falling speed of the sphere. The 
Reynolds number based on sphere size is of order 10*^. 
Thus the sphere sheds turbulent wakes, and this gives 
rise to the stochastic motion. 



Sphere 


Pe (g/cm=^) 


D (cm) 


king-pong 


0.122 


4.41 


ping-pong 


0.146 


3.80 


wood 


0.987 


1.27-3.70 


polypropylene 


1.14 


0.56-2.54 


nylon 


1.56 


0.63-2.54 



TABLE I: Inertial mass density and diameter for the various 
spheres. The ping-pong and king-pong balls are both hollow 
plastic spheres, with a 0.4 mm shell thickness; all others are 
solid. 



IV. EFFECTIVE TEMPERATURE 

In this section we begin reporting on how the various 
terms in the Langevin Equation scale with system pa- 
rameters. The first is the effective temperature, given by 
the mean-squared speed as kT = ■me{v'^)/2. Data for the 
mean-squared sphere speed is shown as a function of air 
speed u in Fig. ^ for various types of sphere. In all cases, 
the data are inconsistent with the simplest dimensionally 
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FIG. 1: The mean-squared speed of a rolling sphere vs the 
speed of the upflow of air, for several types of sphere as 
labelled. The solid lines are a best fit to cubic behavior, 
(u^) ~ u"^. The data are not consistent with the dimensionally 
simpler scaling (u^) ~ u^ , shown as a dashed line. All data 
are for the full 12-inch sieve, except the small right triangles 
for ping-pong balls in cells of smaller radii. 



correct scaling, (w^) ^ m^. Rather, the mean-squared 
speed appears to scale as the cube of the air speed. Thus 
u^ / {v'^) has units of speed and presumably depends on 
physical characteristics of the sphere, the fluidizing air, 
and gravity. 

To uncover the full scaling we first proceed by 
dimensional analysis. Assuming that (w^) de- 

creases vifith increasing sphere density, the combination 
{Pair I Pe^u^ / {v^) IS thc important characteristic speed, 
where the exponent a of the density ratio is to be de- 
termined. We can conceive of only three possibilities for 
the origin of this characteristic speed: the speed of sound, 
34,000 cm/s; a speed set by gravity and the sphere size, 
\fgD\ and a speed set by air viscosity and sphere size, 
r\lD. To investigate, we compare these possibilities with 
data for the characteristic speed vs sphere size in Fig. |21 
for several integer values of a. We find that the best 
data collapse is attained for a — 2. For that case the 
value and functional form of the characteristic speed are 
both consistent with \fgD. Adjusting the numerical pref- 
actor to best match all the data, we thus find that the 
mean-squared speed of a sphere is given by 



0.7 



Pai 



(4) 



This observed scaling of the mean-squared sphere 
speed is consistent with a simple model of the stochas- 
tic motion of the sphere being driven by turbulence in 
the air. The idea is to balance the rate Pin at which 
kinetic energy is transferred from the air to the sphere 
with the rate Pout at which energy is dissipated by drag. 
Ignoring numerical factors, the latter is the character- 
istic drag force times the characteristic sphere speed: 
Pout — (pe-D^row)u, using the notation of Section^ The 
former is Pin — {peD'^Sv^)^o, where the term in paren- 
thesis is the kinetic energy change due to the shedding 
of a wake and 70 = u/D is the rate at which wakes are 
shed. Assuming that the wake size scales with sphere 
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FIG. 2: Scaling of the characteristics speeds with the sphere 
diameter, D, for several different spheres as labelled. The 
mean-squared speed {11^ ) of the sphere is proportional to the 
cube of the air speed, u^; therefore, the ratio of these quan- 
tities is a characteristic speed that reflects both the sphere 
density and the dissipation mechanism. The data collapse is 
best when u^ /{v^) is multiplied by the square of the density 
ratio. Then the value and form of the characteristic speed are 
both consistent with \/gD, indicating that rolling friction is 
the dominant dissipation mechanism. 



size, momentum conservation gives Pe5v = PairU. The 
numerical prefactor is nontrivial, since it must depend 
on the ratio of wake to sphere size and also on the frac- 
tion of momentum in the plane of the sieve, transverse 
to the average air fiow direction. Combining all these el- 
ements, the balance of power input with power output is 
"^u^/D = Pe'^{v'^)To- This is identical to our data on 



P. 



the mean-squared speed, Eq. Q), provided that the drag 
amplitude scales as Fo ex ^ gjD and that the memory 
decay rate scales as 70 ex ujD. Next we demonstrate 
that these provisos both hold true. 



V. DRAG AND RANDOM FORCES 

Recall from Eqs. (jll,S|l in Section|n]that both the drag 
and random forces are specified by an exponential mem- 
ory kernal, r(t) = Fo7o exp(— 7ot). In Ref. 15] we found 
consistent values for Fo and 70 from two different meth- 
ods. The first was from the velocity autocorrelation func- 
tion using thc Langevin equation. The second was from 
the amplitude and phase of the average response to a 
small sinusoidal rocking of the entire apparatus at vari- 
ous frequencies. Here we employ the former method for 
both the ping-pong and king-pong balls, as a function 
of air speed. The results are shown in Fig. |21 made di- 
mensionless according to the expectations of Section ITVI 
Specifically, the top plot demonstrates that the drag am- 
plitude behaves as expected: 



Fo = 0.17 V^. 



(5) 



The importance of g suggests that rolling friction is the 
dominant source of drag, as opposed to shear or compres- 
sion of the air. Perhaps we may identify the numerical 
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FIG. 3: Amplitude To and decay rate 70 of the memory ker- 
nal, r(i) = ro7o exp(— 7oi), as a function of air flow speed, 
for two different spheres as labelled; these quantities are ren- 
dered dimensionless by appropriate factors of sphere diame- 
ter, gravitational acceleration, and air flow speed according to 
expectation. The dashed lines represent average values, 0.17 
in the top plot and 0.11 in the bottom plot. 



prcfactor of g as a coefficient of friction. To — ^ ^g/ D 
with jj, — 0.03. It would be interesting to investigate how 
/i changes with mesh size and ball roughness. The bot- 
tom plot of Fig. 121 demonstrates that the memory decay 
rate also behaves as expected: 



7o 



Q.llu/D. 



(6) 



This is consistent with earlier visualization and pressure 
fluctuation studies, which found that the vortex shedding 
frequency is 0.15m/D for Reynolds number in the range 
10^ - 10^ [HiOl. Here, Eq. © means that a new wake 
is shed every time the air flows a distance of about nine 
sphere diameters; equivalently, the Strouhal number is 

StEE7oD/u = 0.11. 

We emphasize that while the results of Eqs. (|5I6|I di- 
rectly specify the drag force, they also specify the random 
force via the Fluctuation-Dissipation Relation Eq. ^. 
The random driving and the drag forces are different as- 
pects of the same physical interaction between the sphere 
and the turbulence it generates in the air. To recap, the 
random force has Gaussian components and an exponen- 
tial temporal autocorrelation, 

(F,(t') • ¥r{t)) = 2mekTToio exp[7o(t - t')], (7) 

where kT = TOe(w^)/2 is specified by Eq. Q). 

VI. BALL- WALL INTERACTION 

The potential V{r) is the only part of the Langevin 
Equation not yet discussed. This can be deduced from 
the radial position probability function, P{r), using prin- 
ciples of statistical mechanics. Namely, the probability 
to find the sphere in a thin ring of radius r is proportional 
to the ring radius times a Boltzmann factor, 
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FIG. 4: Interaction potential between a variety of spheres 
and the walls of the container, vs distance from the center of 
the cell scaled by cell radius. The top data set is taken at 
constant air flow speed, while all others are taken at con- 
stant (full) cell size. The dashed curves represent a har- 
monic potential, Vh{r)/kT = 30{r / Rceii)^ ■ The solid curves 
represent an empirical average of all the data, V{r)/kT — 
30{r/Rceii)'^/[i- + 2{r/Rceii)'^]- An arbitrary constant offset 
was added in order to separate the different datasets. 



where kT is the effective temperature discussed in Sec- 
tion IIVI In Ref . [13 , the sphere was found most fre- 
quently near the center of the cell such that the x 
and y distributions were nearly Gaussian and P{r) « 
(2r/(r^)) exp(— r^/(r^)). This means that the interac- 
tion potential is nearly harmonic, V{r) « Kr^/2. The 
value of (r^) gave a spring constant that was verified by 
an auxiliary mechanical tilting measurement. Here, we 
examine the shape of the potential more closely, and we 
explore its behavior as a function of system parameters. 
Radial position probability data for all runs are con- 
verted to the interaction potential via Eq. ©, and dis- 
played altogether in Fig. 0| The potential is left in units 
of kT, and the radial position is scaled by the cell radius, 
Rceih for clarity. Remarkably, the potential is given by 
the same empirical form independent of sphere size, cell 
radius, and air flow speed: 



V{r)/kT = 



1 + {r/RceiiY- 



(9) 



P{r) cxrexp[-y(r)/A:T], 



(8) 



In particular, the rms radial position of a sphere is always 
set by the ceh size, a/(^ = (0.20 ± 0.01)Rceii- The 
harmonic form of the potential also softens away from 
the center. It actually becomes attractive very close to 
the walls, strong enough to occasionally trap an unwary 
sphere that wanders too far from home. 

The geometric scaling of the potential with cell size, 
independent of air flow speed, leads us to believe that 



the origin of the behavior lies in the interaction of the 
shed vortices with the boundary of the cell. This is bol- 
stered by other observations as well. First, even very 
slight imperfections in the circularity of the cell can break 
the radial symmetry of the position distributions. Sec- 
ond, placing a hand or other object downstream from the 
sphere affects its position distribution as well. Evidently, 
the vortex street is connected to the sphere such that 
force can be exerted on the sphere via perturbation to 
the vortices. 

One possible picture for how the vortex street senses 
the wall is that the transverse extent of the vortices 
grows linearly with distance downstream. Then the 
sphere could sense its position from the height at which 
the expanding vortices hit the boundary. See Figs. 31, 
55, 56, 172, and 173 of Ref. ^ for photographs of 
the vortices behind various objects at approximately the 
same Reynolds number as here. Another possible pic- 
ture is that the background flow, while homogeneous 
near the sieve, develops large-scale structure downstream 
that grows from the edge inwards. Then the sphere 
could senses its position from from the height at which 
its vortices merge with the dome of turbulent structure 
above. See Fig. 152-153 of Ref. [2l| for photographs of 
the isotropic turbulence behind a grid and it's evolution 
downstream. 

We performed a few tests in attempt to clarify the 
physical pictures. First we increased and decreased the 
wall height to considerable extent. This had no influ- 
ence on the sphere position statistics, which seems to 
rule out the growing-vortex scenario. Our second test 
was to stretch a fine netting across the top of the sieve. 
We hoped that this would affect the rate of vortex shed- 
ding or the way the vortex street is connected back to 
the sphere. However, it had no influence on the sphere 
position statistics either. Thus, we must leave the ori- 
gin of the geometric nature of the sphere- wall interaction 
potential as something of a mystery. Flow visualization 
may be helpful. We close by emphasizing that, whatever 
its origin, the sphere is repelled by the cell wall in a way 
that, remarkably, can be described by a potential energy 
and a corresponding conservative force. 



VII. BALL-BALL INTERACTION 

In the remainder of this paper we report on the air- 
mediated interaction between two spheres rolling in the 
same nearly-levitating upflow of air. Throughout, we set 
the air flow to 280 cm/s, before adding spheres. As above, 
we shall see that this may be studied using position prob- 
ability data and statistical mechanics. And just as for 
the ball-wall force, we shall see that the ball-ball force is 
repulsive. Naively one might expect a Bernoulli-like at- 
traction, just as when air is blown between two objects. 
However it's immediately obvious from visual inspection 
that here the two spheres repel. Only rarely do they col- 
lide, with physical contact between their surfaces; they 
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FIG. 5: Speed and radial position probability functions for 
one and two spheres rolling in the same upflow of air. Note 
that the speed distributions are Gaussian in both cases, as 
shown by the dashed curves. By contrast, the radial posi- 
tion function becomes non- Gaussian when a second sphere is 
added to the system. 



never stick; they accelerate apart after close approach. 

To begin we display speed and radial position proba- 
bilities in Fig. 13 The light curves are for a single ball 
in the same air flow, for comparison. As above and in 
Ref. [13, the x and y components of velocity and posi- 
tion are all Gaussian. When a second sphere is added, we 
verify that the velocity and position distributions remain 
radially symmetric and identical for each sphere. The 
top plot of Fig. [SI demonstrates that the average speed 
distribution of the two spheres remains Gaussian. Thus 
the mean-squared speed can be used to define an effective 
temperature, as before for a single sphere. However, this 
temperature increases when a second sphere is added, 
even though the flux of air remains unchanged. Evidently 
Eq. Q holds in detail only for a one-sphere system. The 
reason may be that, due to a decrease in free area, the air 
flow speed around the two spheres is greater than when 
only one is present. It may also be that the process of 
energy injection via vortex shedding is altered. The bot- 
tom plot in Fig. |31 demonstrates that the radial position 
probability becomes non-Gaussian when a second sphere 
is added. Each sphere spends less time in the very center 
of the cell, due to mutual repulsion, with the rms radial 
position increasing from 2.8 cm to 4.8 cm when a second 
sphere is added. As before, the spheres still are repelled 
from the cell wall as though in a harmonic trap. 

For statistical mechanics to be useful for studying the 
sphere-sphere repulsion, it is required that the velocity 
components be Gaussian as demonstrated above. It is 
also required that there be no correlation between the in- 
stantaneous velocities of the two spheres. To check this, 
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FIG. 6: Velocity correlations between spheres i and j, for 
one- and two-sphere systems as labelled. Note that the cross- 
correlation vanishes at t = 0, which is required if statistical 
mechanics is to be invoked. 



we compute temporal velocity correlation functions and 
plot the results in Fig. El The velocity autocorrelation 
for a single sphere alone in the cell is shown by a light 
curve, for comparison. The velocity autocorrelation for 
each sphere, when two are present, is shown by a heavier 
curve. It decays over the same time scale as the one- 
sphere autocorrelation, though the oscillations are less 
pronounced. The cross correlation between the velocities 
of the two spheres is shown by a dashed curve. It too 
oscillates and decays over the same time scale as the au- 
tocorrelations. But, crucially for us, it vanishes at r = 0. 
Thus the instantaneous equal-time velocities of the two 
spheres are indeed uncorrelated as required. 

With the above preliminaries established, we may now 
exploit the principles of statistical mechanics in order 
to deduce the sphere-sphere interaction potential Vssip), 
where p is the distance between the centers of the two 
spheres. The idea is to compute the sphere-sphere sepa- 
ration probability in terms of both the overall harmonic 
confining potential and the unknown Vss{p)- This is ac- 
complished by summing the Boltzmann factors for all the 
ways of arranging the spheres with the desired separation: 



P{p) (x / dxdydO exp 



K{ x^ +y^ 



(10) 



{x + pcos9f + {y + p sin 9f j /kT 
xexp[-Vss{p)/kT] 



oc exp 



-[lKp^ + VUp)]/kT 



(11) 



One may differentiate this expression to show that the 
peak in P{p) is where —dVss/dp = Kp/2, which is a 
statement of force balance when each sphere is p/2 from 
the center of the cell. Since the spring constant K is 
known from the one-sphere experiment, and since the 
temperature kT is known from the mean-squared speeds, 
the functions P{p) and Vss{p) may be deduced one from 
the other. 

The separation probability P{p) is readily found from 
the video data for the position of each sphere vs time. 
Results are displayed by a dashed curve on the right axis 
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FIG. 7: Interaction between two spheres as a function of their 
separation p. Using Eq. Ijll^ . the potential Vss{p) is inferred 
from the separation probability P{p); these functions are both 
shown in the upper plot. Systematic uncertainty in Vss{p) is 
indicated by dotted curves; it is due to statistical uncertainty 
in the value of K. The repulsive force is shown in the bottom 
plot, as obtained both from Ks(p) and from an auxilliary 
mechanical measurement. 



of the upper plot in Fig.d The probability rises abruptly 
from zero at a separation equal to the sphere diameter. 
It reaches a peak near p — 7 cm, and then gradually 
decays again toward zero. The sphere-sphere potential 
Vss{p) can then be obtained from P{p) using Eq. Hll|l . 
Results are shown by a solid curve on the left axis of the 
upper plot in Fig.[3 The precipitous drop of Vss{p) near 
contact indicates a hardcore repulsion. The more gradual 
drop at larger separations indicates a softer repulsion. 

The actual force of repulsion may be found by differ- 
entiating, Fss{p) — —dVss/dp. Results are shown by 
the solid curve in the lower plot of Fig. There is a 
hardcore repulsion, followed by a nearly constant-force 
repulsion when the sphere centers are separated by more 
than two diameters. Expressing the interaction in terms 
of a force allows us to perform a check using an auxiliary 
mechanical measurement of the response to tilting the 
entire apparatus by a fixed angle 9 away from horizontal. 
This causes a constant component of gravity, mgsm9, 
within the plane and breaks the radial symmetry; note 
that here m is the true mass, not the effective inertial 
mass. Then we measure the probability (t>{x, y) for find- 
ing a sphere at a given position, where the origin of the 
coordinate system is at the center of the cell and where 
gravity acts in the +y direction. This probability has 
two peaks, at coordinates (±p/2,yp), separated by dis- 
tance p. Assuming only that the wall repulsion acts in 
the radial direction, the statement of force balance at the 
peaks of 4>{x, y) gives the sphere-sphere repulsive force as 



*(P) = ( ^P/yp ) "^ffsin6'. 



(12) 



In practice, to achieve a wider range in separations, we 
tilt the apparatus by 0.013 rad and use ceUs of three dif- 
ferent diameters: 20, 25, and 30 cm. Observations then 
give the repulsive force at three different separations as 
shown in the lower plot of Fig. [7| Evidently the agree- 
ment with the results from statistical mechanics is very 
good. This gives confidence in the use of statistical me- 
chanics to deduce the full form of the repulsive sphere- 
sphere interaction. 



VIII. CONCLUSION 

We have exploited the thermal-like behavior of a single 
gas-fluidized sphere to deduce the nature of the forces 
dictating its motion. All these forces are mediated by 
turbulence in the gas, but can be decomposed into dis- 
tinct contributions. Due to randomness in the shedding 
of turbulent wakes, there is a rapidly varying random 
force specified by Eas. l|5l7|l . By virtue of the Fluctuation- 
Dissipation Relation, and Eas. Hll2|l . these results also 



fully specify a velocity-dependent drag force that damps 
rolling motion. The apparent interaction of the wakes 
with the cell boundary gives rise to a nearly-harmonic 
force that keeps the rms sphere position at about one 
fifth the cell radius, no matter how the system parame- 
ters are changed. The effective temperature, set by the 
mean-squared speed in Eq. Q), is a key parameter in 
these forces. When a second sphere is added, the ther- 
mal analogy still holds and these forces change only in 
detail. In addition, there is a gas-mediated repulsion act- 
ing between the spheres that is nearly constant beyond a 
few diameters of separation and that grows stronger near 
contact. 
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